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Abstract. In the framework of the relativistic mean- field theory, we have considered the equa- 
tion of state of superdense nuclear matter, taking into account an effective scalar-isovector 
(5-meson field. The effect of the i5-meson field on the characteristics of a Maxwell-type quark 
phase transition has been studied. The quark phase is described with the aid of the improved 
version of the MIT (Massachusetts Institute of Technology) bag model, in which interactions 
between the u, d, s quarks inside the bag are taken into account in the one-gluon exchange 
approximation. For different values of the bag parameter i?, series of neutron star models with 
a quark core have been built. Stability problems for neutron stars with an infinitesimal quark 
core are discussed. An estimate is obtained for the amount of energy released in a catastrophic 
transformation of a critical neutron star to a star with a finite-size quark core. 
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1 Introduction 

The properties of such compact objects as neutron stars (NS) depend functionally on the equation 
of state (EOS) of matter in a sufficiently widerange of densities. A knowledge of the structure 
characteristics and the constituent composition of matter at extremely large densities is a necessary 
condition for revealing the physical nature of the NS internal structure and integral parameters. 
One of the effectively used theories, sufficiently adequately describing the properties of nuclear 
matter as a system of strongly interacting baryons and mesons, is the relativistic mean field (RMF) 
theory [1, 2]. This theory has made possible to obtain results which satisfactorily describe the 
structure of finite nuclei [3], the EOS of nuclear matter [4], and the features of heavy ion scattering 
[5]. Inclusion of a scalar-isovector 6 meson into the scheme and a study of its role for asymmetric 
nuclear matter in the low density range has been conducted in [6, 7, 8]. An objective of this paper 
is to study the EOS of superdense nuclear matter in the framework of the RMF theory and an 
investigation of changes in the parameters of the first-order phase transition due to inclusion of 6 
meson exchange. Using the EOS obtained, we calculate the integral and structure characteristics 
of NS with a quark core. 
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2 Properties of Neutron Star Matter 



2.1 Nuclear Matter 



The nonlinear Lagrangian density of an interacting multiparticle system consisting of nucleons 
and an isoscalar-scalar a meson, an isovector-scalar 6 meson, an isoscalar-vector u; meson and an 
isovector- vector p meson as well as free electrons has the following form in quantum hadrodynamics 

(QhdS 



C = 'lpN 



+ 



\ {d^cr{x)d''a{x) - m^a{xf^ - ^m^ {g^cr{x)f - ^ ig^cr{x))^ 



+ 



(1) 



Here, g^, g^j, gs and gp denote the coupling constants of the nucleon with the corresponding meson. 
In the RMF theory, the meson fields a (x), oj^ (x), 5 (x) and (x) are replaced by the (effective) 
fields a, (D^, 5, p^. Re-denoting the meson fields and coupling constants according to 



9a<y = o- , g^ujQ = a; , gs 



(3) 



(2) 



{gcj/m„f = , {gu^/m^f = , {gs/msf = as , {gp/nipf = ap 
and introducing the asymmetry parameter 

a= {un- np)/n , 

one can present the field equations in the form 

V {risp {n, a) +nsn {n, a) - brriNa'^ - ca^^ , 



a = an 



uj = au,n . 



5 = as {usp (n, a) -Usn {n, a) ) , 



p = --apua 



where 



nsp{n,a) = ^ J 



ml (a, 5) r, 

' k^dk , 



^k'^ + m* (fj, 6) 



We will use the natural system of units h — c = 1 



(3) 
(4) 

(5) 

(6) 
(7) 

(8) 
(9) 
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Table 1: Constants of the theory without {(JLop) and with (aujpS 





Oct, fm^ 


auj, fm^ 


as, fm^ 


Op, fm^ 


b, fm ^ 


c 


atop 


9.154 


4.828 





4.794 


1.654 • 10-2 


1.319 • 10-2 


aiop6 


9.154 


4.828 


2.5 


13.621 


1.654 • 10-2 


1.319 • 10-2 



the S meson field 



.2,_ \ 1/3 



fcF(n)=(^l . (11) 

The effective masses of the proton and neutron are determined by the expressions 

m*{a,6) = mN— a — 6, m*^{a,5) = rriN — cf + 6. (12) 

If the constants a^^ and Up are known, equations (6) and (8) determine the functions lo (n) and 
p{n,a). Moreover, a knowledge of the other constants a^, as, b, and c makes it possible to solve 
the set of equations (5), (7), (9), (10) in a sclf-consistcnt way and to determine the remaining 
two meson field functions a (n, a) and 5 {n, a). The standard QHD procedure makes it possible to 
obtain expressions for the energy density e (n, a) and pressure P (n, a) [9]. 

The constants of the theory arc numerically defined in such a way as to reproduce the em- 
pirically known values of such nuclear characteristics at saturation as the bare nucleon mass 
rriN = 938.93 McV, the parameter 7 = m*j^/mM = 0.78, the saturated nuclear matter baryon 
number density no = 0.153 fm~^ , the binding energy per baryon /o = —16.3 MeV, the compress- 
ibility modulus K = 300 MeV, and the asymmetry energy Eiyln = 32.5 MeV. To reveal the role 
of taking into account the S meson, we have used the value as = 2.5 fm2 (see [10]). Table 1 shows 
the values of constants of the theory obtained without taking into account the interaction channel 
due to the isovector scalar S meson {aup) and taking it into account {aupS). 

Fig. 1 represents a three-dimensional picture of the energy per baryon Ei,{n,a) = £NM/n 
vs. the baryon number density n and the asymmetry parameter a in the case of a /3-equilibrium 
charged "npe"-plasma. The lines correspond to different fixed values of charge per baryon. The 
thick line corresponds to a /3-equilibrium, electrically neutral "npe" -matter. The lower surface 
corresponds to the "acop" model, the upper one to the "aiop6" model. It is seen that inclusion of 
the 5 meson field increases the energy per nucleon value, and this change is enhanced as the nuclear 
matter asymmetry parameter increases. The results of our analysis show that the scalar-isovector 
(5-field leads to an increased stiffness of the nuclear matter EOS due to splitting the effective masses 
of the proton and the neutron thus increasing the asymmetry energy. The EOS obtained by us in 
the normal nuclear density range is matched to the well known Baym-Bethe-Pethick (BBP) EOS 
[11]. 
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Figure 1: A three-dimensional representation of the energy per baryon as a function of the 
baryon number density n and the asymmetry parameter a in the case of a /3-equihbrium charged 
"npe"-plasma. The upper surface corresponds to the ^^aujp5" model, the lower one to the ^^awp" 
model. 

2.2 Strange Quark Matter 

To describe the quark phase, we have used the improved version of the MIT (Massachusetts 
Institute of Technology) bag model [12], in which interactions between the u, d, s quarks inside the 
bag are taken into account in the one-gluon exchange approximation [13]. The quark phase consists 
of the three quark flavors u, d, s and electrons in equilibrium with respect to weak interactions. 
For the quark masses we have used the values m„ = 5 MeV, = 7 MeV and = 150 MeV, 
while for the strong interaction constant we have taken the value as = 0.5. The calculations were 
conducted for different values of the bag parameter B in the range of 60 to 120 MeV/fm^. 

2.3 Phase Transition to Quark Matter under Constant Pressure 

We assume in this paper that the nuclear to quark matter conversion is the ordinary first-order 
phase transition described by the Maxwell construction. A separate paper will be devoted to a 
study of alterations in the characteristics of a phase transition with mixed phase formation [14] 
due to inclusion of the 5-meson field as well as the effect of these alterations on the integral and 
structure parameters of hybrid stars. In the case of an ordinary first-order phase transition, it is 
supposed that both the nuclear and quark matter, taken separately, is electrically neutral, and at 
a certain value of the pressure, Pq, corresponding to co-existence of the two phases, the baryonic 
chemical potentials of the two phases coincide. Table 2 represents the parameter values of such a 
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Tabic 2: Parameters of the Maxwell phase transition for different values of the "bag" constant B 
without taking into account the 6 meson field {aujp)and with it {aujp5). 



B 

MeV/fm^ 




fJ-b 
MeV 


UN 




Po 

MeV/fm^ 


En 

MeV/fm^ 


MeV/fm^ 


A 


60 


atop 


965.4 


0.1220 


0.2826 


1.965 


115.8 


270.9 


2.299 


60 


aujpS 


965.9 


0.1207 


0.2831 


2.11 


114.5 


271.4 


2.327 


69.3 


aojp 


1037 


0.246 


0.3557 


15.57 


239.6 


353.4 


1.385 


69.3 


acopS 


1032 


0.2241 


0.3504 


13.84 


217.5 


347.9 


1.504 


80 


aup 


1142 


0.3792 


0.4819 


48.54 


384.5 


501.8 


1.159 


80 


aup6 


1119 


0.3276 


0.4525 


37.95 


328.8 


468.6 


1.278 


100 


atop 


1298 


0.5506 


0.7175 


121.3 


593.4 


810 


1.133 


100 


acopS 


1257 


0.4746 


0.6497 


93.30 


503.3 


723.5 


1.213 


120 


atop 


1396 


0.6512 


0.8975 


180.2 


728.9 


1073 


1.18 


120 


aup5 


1354 


0.5729 


0.8165 


143.9 


631.7 


961.4 


1.24 



phase transition for five different values of the "bag" parameter B both with and without the S 
field. In this table, pb is the baryonic chemical potential at the phase transition point, n^y and uq 
are the baryon number densities of the nuclear and quark matter, respectively, at the transition 
point, En and eq are the energy densities, Pq is the phase transition pressure, and A = eq/(en+Po) 
is the density jump parameter. It can be seen that inclusion of the S interaction channel results in 
decreased values of the transition pressure Pq and baryon number densities of both phases at the 
phase transition point, njv and uq. Meanwhile, the value of the jump parameter A increases. 

3 Models of Neutron Stars With a Quark Core 

Using the EOS of nuclear matter obtained in the previous section we have integrated the Tolman- 
Oppenheimer-Volkoff (TOV) set of equations [15] and obtained both the structure functions e(r), 
P{r) and m{r) and the stellar integral parameters: the gravitational mass M and the radius R for 
different values of the central pressure Pc- Fig. 2 shows the M{R) dependence for different values of 
the bag parameter B. It is seen that the purely nucleon EOS leads to a maximum mass of ~ 2.2 Mq. 
The advent of a quark phase in the NS diminishes this value up to ~ 1.75 Mq. According to [16], in 
the case A > Xcr = 3/2, an infinitesimal core of the new phase is unstable. Our analysis shows that, 
for B < 69.3 MeV/fm^, the density jump parameter has the values A > Xcr, and NS configurations 
with infinitesimal quark cores will be unstable. In the latter case, there is a nonzero minimum 
value of the quark core radius of a neutron star. Accretion of matter on a critical neutron star will 
then result in a catastrophic (jumplike) transformation of the star, forming a star with a quark 
core of finite size. The catastrophic transformation with finite quark core formation at the stellar 
center will release an enormous amount of energy, compared with energy release at a Supernova 
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Figure 2: The stellar mass M vs. radius R for different values of the bag parameter B in the ^^acopS" 
model. The values of B in MeV/fm^ are shown near the critical configurations corresponding to 
quark phase formation. 

explosion. Our calculations have allowed an evaluation of this energy. Thus, for B = 60 McV/fni^, a 
star of mass M ~ 0.24 Mq and radius R ~ 16.77 km forms, in a jump-like manner, a star of radius 
R 13.95 km and a quark core, having a core mass Mcore ~ 0.087 Mq and radius Rcore ~ 4.38 km. 
The catastrophic restructuring process releases an energy Econv ~ 4 10^^ J. For B > 90 MeV/fm^, 
the quark degrees of freedom make the neutron star unstable, and configurations with a quark core 
are absent. In the case 69.3 < S < 90 MeV/fm^, configurations with arbitrarily small quark cores 
are stable. 

4 Conclusion 

In the framework of the RMF theory, we have considered the EOS of NS matter and obtained that 
inclusion of the isovcctor-scalar (5-mcson field leads to nonncgligible changes in the characteristics of 
quark first-order phase transitions. For different values of the bag parameter B, we have constructed 
neutron star models with quark cores. The results have shown that, depending on value of this 
parameter realized in the Nature, the following criteria may take place. In the case B > 90 
MeV/fm^, the quark phase leads to a total neutron star instability, and stars with quark cores are 
absent. In the case 69.3 < B < 90 MeV/fm^, stable NS with arbitrarily small quark cores may 
exist. If B < 69.3 MeV/fm^, then configurations with small quark cores are unstable. In the latter 
case, accretion of matter will result in a catastrophic NS transformation after which a quark core 
of finite size will be formed. 
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